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Generalized parton distributions and double distributions for qq pions 
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We consider two simple covariant models for pions (one with scalar and the other with spin-i 
constituents). Pion generalized parton distributions are derived by integration over the light-cone 
energy. The model distributions are consistent with all known properties of generalized parton 
distributions, including positivity. We also construct the corresponding double distributions by 
appealing to Lorentz invariance. These ostensibly constructed double distributions lead to different 
generalized parton distributions that need not respect the positivity constraints. This inconsistency 
arises from the ambiguity inherent in defining double distributions in the standard one-component 
formalism (even in the absence of the Polyakov- Weiss term). We demonstrate that the correct model 
double distributions can be calculated from non-diagonal matrix elements of twist- two operators. 

PACS numbers: 13.40.Gp, 13.60.Fz, 14.40.Aq 



I. INTRODUCTION 

In recent years generalized parton distributions 
(GPDs) have generated a considerable amount of at- 
tention. These distributions stem from hadronic ma- 
trix elements that are both non-diagonal with respect 
to hadron states and involve quark and gluon operators 
separated by a light-like distance. Thus physics of both 
inclusive (parton distributions, e.g.) and exclusive (form 
factors, e.g.) reactions is contained in the GPDs. At the 
leading- twist level, these new structure functions describe 
the soft physics of a variety of hard exclusive processes 
(see the reviews |^). 

Since light-like correlation functions are involved in 
the description of deeply virtual Compton scattering 
(DVCS), e.g., there exists a simple decomposition of these 
matrix elements in terms of the light-cone Fock space 
wave functions of the initial and final states ^. This 
representation of GPDs is ideal for physical intuition, 
however, comparatively little has been done to show that 
the light-cone wave function representation is consistent 
with the reduction properties required of the generalized 
parton distributions. Below we undertake a simpler task 
of presenting covariant models for the pion which respect 
the properties of GPDs. Albeit simple, these models il- 
lustrate the utility and physicality of the light-cone Fock 
representation as well as provide a guide to understand- 
ing how the reduction relations arise in this formalism 
which will be useful when non-perturbative solutions for 
the Fock components in QCD become available. The 
scalar constituent model which we consider is merely the 
triangle graph with point-like vertices (this has been con- 
sidered in see also 0). The spin-i model is based on 
an earlier exactly soluble, (1 -I- l)-dimensional light-front 
model 0. Extension of this model to (3-1- 1) dimensions, 
which involves regularizing the divergent light-front cur- 
rent, has been done in 

Another approach is to use the formalism of double 
distributions (DDs) |lO| . This formalism elegantly ex- 
plains the polynomiality conditions required of GPDs and 
thus gives one the ability to construct models consistent 



with known properties — although insight into model con- 
struction has often been limited to factorization Ansdtze. 
Recently two-body light-front wave function models of 
the pion have been used to obtain GPDs jl^ based on 
DDs. Without modifying the quark distribution, the 
resulting GPDs violate the positivity constraints [ p^ . 
This inconsistency is attributed to missing contributions 
from non-wave function vertex diagrams which are absent 
when one uses non-covariant vertices. In general these di- 
agrams are a substitute for higher Fock components, see 
]l3[ . The covariant models used here, however, allow us 
to test the uniqueness of this ostensible construction. 

Indeed we find that appealing to Lorentz invariance is 
only enough to determine one component of the double 
distribution in the two-component formalism (even for C- 
odd distributions, where the Polyakov- Weiss D-term 
is absent) . We show that in the scalar constituent model, 
missing the second component leads to incorrect GPDs. 
The same is true for the spin-i constituent model, where 
additionally the positivity constraint is violated. The 
correct DDs can be calculated from non-diagonal matrix 
elements of twist-two operators which we demonstrate 
for both models. On the other hand, exploiting the am- 
biguity inherent in defining DDs (which is akin to gauge 
freedom p^) one can generate infinitely many different 
GPD models which share the same form factor and quark 
distribution as well as satisfy polynomiality (and likely 
positivity) . 

The org anization of the paper is as follows. First in 
section ||, we explicitly derive the GPD for the scalar 
triangle diagram with point-like vertices. Next we show 
the DD for this model extracted from the form factor 
in the Drell-Yan frame leads to incorrect GPDs. Having 
encountered this problem, we rev iew definitions of the 
double distributions in section III. Here we also calcu- 



late the correct DD for the scalar triangle diagram. In 



[V, we present the model for the spin--^ case. Next in 



section 



IV B, we regularize the current and then extract 
this model's GPD. Although not manifest, this model 
satisfi es polynomiality, which is demonstrated in section 
IV C. Using as a guide, we construct the DD for this 
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FIG. 1: Covariant triangle diagram for the pion electromag- 
netic form factor 



model in section IVD. Similar to section O, this one- 



component DD too gives rise to a different GPD than 
the hght-front projection. Additionally positivity is not 
satisfied by this one-component DD (section IV E). Wc 
calculate the missing component of the DD from matrix 
elements of twist-two operators in section [[VF|. Lastly 
we conclude with a brief summary (section] 



II. PION WITH SCALAR CONSTITUENTS 

For the pion model with scalar constituents, we choose 
the point-like Bethe-Salpeter vertex T{k,P) = 1, where 
the couphng constant is assumed to be absorbed into the 
overall normalization. Furthermore, we choose derivative 
coupling of the photon to charged scalar particles. 

The pion electromagnetic form factor for this model 
can be calculated from the Feynman triangle diagram. In 
order to derive the GPD, however, we need to choose the 
kinematics specified in Figure |^ with k as the momentum 
of the struck quark. Using the stated pion vertex and 
taking the plus-component of the current^, we have 



< [(fc-f A)2-m2+ie]"^[(P-fc)2-rn2+ie]"\ (1) 



where the momentum transfer ist = and the skewness 
is defined relative to the initial state A+ = — CP^ < 0. 
Physically ( plays the role of Bjorken variable for DVCS. 
Additionally we work in the frame where = 0. 

To turn Eq. (|l|) into an expression for the GPD 
H(x, C>^)j insert 6(k^ /P^ — x) to fix the momentum 
of the struck quark and keep C 7^ 0. This forces 



Fit) 



H[x, t)dx. 



(2) 



Lastly we integrate over k to project onto the light cone. 
Doing the contour integration to extract H{x,C,t) in 



For any vector , we define the light-cone variables = 
-I=(a0±a3). 



Eq. (|i|), we are confronted with the poles 

'^a — %n X 

k- = p- + {k-p)-,~^^ 

.fcc" = -A- + (fc + A)-„-i^ 



(3) 



^ I 2 

where the on-shell energies are p^^ = ^ 2^+™ ^'^^ 



abbreviation 



is used. Thus the non- vanishing 



contribution to the integral is 

27ri 6'[a;(l-a;)](6i(a;-C)R-cs(fc-)-6l(C-2;)Res(fc-)), (4) 

which leads to 

{I - C/2)H{xX,t) = e{x - (:)H,{xX,t) 

+ eiC-x)H2{x,(:,t). (5) 

Using k'^ — k-'- +{1 — x') A.^ for the relative transverse 
momentum of the final state, the functional forms are 



Hi{x,C,t) = {2x-C)\Nf J dk^ Dw(x,k^|M2) 



X Dw(x',k'-L \M^)/x{l~x)x' 



H2{xX,t) - {2x-C)\N\^ I dk^ DwixM 

X Dw(x",k"^|i)/Cx"(l - x")(l - x), 



(6) 



(7) 



where x" = x/(^ and k"^ = k-' 



A are the rela- 



tive momenta of the photon. Additionally, we use the 
replacement 



Dw(a;,k-^ \M^)-^ = 



x{l — x) 



(8) 



which is the propagator of the Weinberg equation |g2| . 

Comments about the GPD H{x, C, t) in Eq. (^ are in 
order. Firstly, the model is covariant and thus the sum 
rule and polynomiality conditions are met (see section 
IVC below for clarification). We have checked this ex- 
plicitly and suitable discussion can be found in |j. 
Secondly, Hi{x,(^,t) appearing in Eq. (^) satisfies the 
relevant positivity constraint for a compound scalar of 
scalar constituents (which appears in the Appendix of 
Jl2t ) which is clear from inspection. 

Consideration of this model was first done from the 
perspective of DDs, see e.g. the toy model of ||]. This 
DD model was revisited recently with derivative coupling 
at the photon vertex in the Appendix of and the same 
DD also appears in To derive the DD for this simple 
model, we appeal to Lorentz invariance, recalling along 
the way the relevant properties of DDs. 

First consider the form factor. In the C = frame, 
Eqs. (^|-0) reduce to the Drell-Yan formula ||2^ via the 
definition in Eq. 



F{t) = \N\'J i^^Dw(x,k^|A42) 



!(l-a;) 

xDw(x,k^+(l-x)A^|M2), (9) 
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FIG. 2: Comparison of covariant GPDs for the scalar triangle 
diagram. The GPDs Eq. (||) (denoted LC) and Eq. (|l|) (DD- 
based) are plotted as a function of x for fixed — 0.7 and 
t — —m? for the mass = m < 2m. We also plot the 
difference between the two curves (5). The area under the 
curves is identically F(-m^) for LC and DD-based GPDs, 
and hence zero for their difference 5. 



where t — — A^'^ . The form factor is Lorentz invariant 
and has a decomposition in terms of the Lorentz invariant 
DD F{x,y;t), namely 



Fit) 



F{x, y; t)dydx. 



The DD satisfies the following relations jlC 
property 

F{x,y;t) (xe{l-x- y), 



(10) 
support 

(11) 



reduction to the quark distribution at zero momentum 
transfer 



pl — X 

q{x) = / F{x,y;0)dy 
Jo 



and is Miinchen symmetric 26 



(12) 



(13) 



Using Eq. (| 
F{x,y;t) 



Fix, y; t) = F{x, l-x-y;t). 

can write F{t) in the form ( p^ ) with 
x\N\^e{l- x-y) 



we 



a;(l - x)M^ -y{l-x- y)t' 



(14) 



which satisfies Eqs. (llj-13|). The ingenuity of DDs comes 
about when we construct the GPD via 

H{xX,t)^ ( f F{z,y;t)6{x-z-Cy)dydz. (15) 
Jq Jo 

In this form the sum rule and polynomiality conditions 
follow trivially. 

In Figure ||, we plot the GPD Eq. (|) as well as the 
GPD derived from DD via Eq. (|l^) . Surprisingly the two 
are different despite the fact both models are covariant 
and posses the same form factor and quark distribution. 
Additionally we plot their difference (S) as a function of 

X. 



III. DEFINING DOUBLE DISTRIBUTIONS 



In this section we define the DDs via their moments following the two-component formalism of p4, |lql . Focusing 
on the ambiguities inherent by defining DDs in this fashion, we will understand why the DD in Eq. (14|) leads to an 
incorrect GPD. Moreover, we shall calculate the correct DD for the scalar triangle diagram from matrix elements of 
twist-two operators. We remark in passing that the two components of the DD (below F and G, or F and Z?-terni in 
the standard formalism) can be viewed as projections of a single function of two variables p6[. 



A. Definitions 

The non-diagonal matrix elements of twist-two operators are more conveniently expressed in variables symmetric 
with respect to initial and final states. To this end we define the average momentum f = (P -f P')'^/2. Let 

D'^ — i(9'' — 9^). Then for a pion of spin-i constituents we have 



( P +A/2\i^{0)jl''iDi'' ■ ■ ■ iL>^"J V^(0)| P -A/2 ) 



2P 



k=Q 



k\{n-k)\ 



-A 



n 



k\{n-k)\ 



(-1)"'' 

A \ Mn-fc+l 

^ 2". 



where the action of ^ • • • ' on Lorentz indices produces only the symmetric traceless part and T-invariance restricts k 
in the first sum to be even and odd in the second. For a pion of scalar constituents, replace 7^ with 2iD^. As it 
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stands there is considerable freedom in the above de- 
composition, e.g. one could rewrite the above with 
k{n— k)Bn,k-i{t) as a contribution to Ank{t)- Carrying 
this out for all fc, puts the bulk in the first term and ren- 
ders the second term proportional only to the symmetric 
traceless part of (n + 1) A's — moments of the Polyakov- 
Weiss _D-term This is the usually encountered form 
of the DD with Z?-term. Calculationally, however, we 
find Eq. (|l6|) is the most useful. 

Summing the moments, we produce the DDs 



Ank(t)^ dp daP^-''a''F{P,a;t) (17) 

/•I /■1-I/3I 

Bnk{t)= dp da l3'^-''a^G{l3,a;t). (18) 

As a consequence of the restrictions on k in the sums, 
the function F(/3, a; t) is even in a while G(/3, a; t) is 
odd. Also for n-even, there is no contribution from the 
ZJ-term to the function G(/3, a; t). 

These functions then appear in the decomposition of 
matrix elements of light-like separated operators 

( P +A/2 I V'(-z-/2) ^ iIj{z-/2) I P -A/2 ) 

= 2Pz dp dae-''^^-'+'°"^-'/^F{P,a;t) 

J-l "'-1+1/31 

fl /-i-l/si 
A z I dp dae-'^^-'+'°"^-'/^G{P,a;t), 

'-1 + I/3I 

(19) 



-1 



where = 0. 

Denoting ^ ~ — A+/2P+, the GPD in symmetric vari- 
ables reads 

Hii.tt)^^^ J dz-e^^*^~ 

X (P+A/2|^^(-z-/2)7+?A(z-/2)|P-A/2). (20) 
Inserting Eq. ( p^ ) into this definition yields 

/•I /-i-l/^l 
H{x,£„t)= dp dad{5i-P-£,a) 

F{P,a;t)+^G{P,a;t)\. (21) 

By integrating over x, we uncover two sum rules: the 
sum rule for the form factor 



(22) 



dp J da F{P,a;t) = F{t) 
and what we call the G-sum rule 

dp I da G{P,a;t) = 0, 



which is trivial since G is an odd function of a. Eq. (|22^ 
has non-trivial consequences however, e.g. it shows the 



(23) 



method employed by |ll| leads only to the F DD in 
Eq. (^. This function integrates to the form factor via 
Eq. (22|) and in the forward limit {^,t 0} reduces to 
the quark distribution (when integrated over a). Thus 
F(P, a; t) should be properly termed the forward-visible 
DD, which encompasses more than just neglecting the D- 
term. From Eq. (|l|), we see that G(/3, a; t) does affect 
higher moments of the GPD. This is the source of the 
discrepancy shown in Figure as we now demonstrate. 



B. Scalar model, revisited 



To derive both F and G DDs for the scalar triangle 
diagram of section O, we must now consider the action 

of the operator 2i^^iZ?^i • • • i_D^"] . This produces a 
factor 



1 



-{2k + A) L^(2/fc + A)^i ■■■{2k + Ay 



(24) 



in the integrand of Eq. (^ . which we now take in the sym- 
metric frame. After the integration over k is performed, 
we are left only with 



1 .1-/3 

dp / daD{P,a;t) 

J -1+13 

(2/3P -aA) [^(/3 P -aA/2y' ■■■{PP -aA/2)^ 



(25) 



where we have used the replacement 

^^'^'"'^^ " m2 - p{l - p)M^ - t[{l - PY - a2]/4' 

(26) 

Using the binomial expansion, we can identify 
F{p,a;t) and G{P,a;t) via Eqs. ^ and (|l|), namely 

F{P,a;t) P 9{l ~ P - a)D{P,a;t)\N\^ (27) 
G{P,a]t) ^ a 6{l - P - a)D{P,a]t)\N\'^ . (28) 

To compare with the results of section ||, we must revert 
to asymmetric variables which is accomplished by {/3 ^ 
x,y (a — p+l)/2}. The denominator common to both 
terms becomes 



m? — x{l — x)M'^ — y{l — x — y)t 



(29) 



Thus we have 



F{x,y;t)=xe{l-x-y)D{x,y-t)\N\^ (30) 

G{x, y- t)^{2y + x~ 1) e{l~x- y)D{x, y; t)\N\\ 

(31) 

Notice the function G{x, y; t) is Miinchen antisymmetric 
which is required because G(/3, a; t) is odd with respect 
to a. 
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To construct the GPD H{x,C,t) we must also convert 
Eq. (^) to asymmetric variables. 

H{x,C,t)= dz dyS{x-z~Cy) 











c 



Fix,y;t) + ^^G{x,y;t) (32) 



We can now plot the GPD in Eq. ( p^) using the two DDs 
in Eqs. ( ^ ) and (^. The result agrees with Eq. (||) 
depicted in Figure^ Moreover, the contribution from 
2Z^G is identically the difference d plotted in the figure. 



IV. PION WITH SPIN-i CONSTITUENTS 

A. Bethe-Salpeter amplitude and pion wave 
function 

For the spin-i model, we choose the trivially qq sym- 
metric Bethe-Salpeter vertex 



(33) 



Here we have assumed only 7^ coupling at the quark- 
pion vertex with coupling constant g, whereas four Dirac 
structures exist This simple coupling is suggested 

by an effective interaction Lagrangian of the form (see, 
e.g. 11) 



£1 = —ig TT . (77^ Tq, 



(34) 



where the coupling constant g ~ ni,/ f^^, with m the con- 
stituent mass and the pion decay constant. Notice 
the (ladder approximation) kernel is independent of light- 
cone time. Thus this model (as well as the scalar triangle 
model in section |l^ are special cases of the instantaneous 
formalism described by jTot in the impulse approxima- 
tion. The relation of the vertex to the Bethe-Salpeter 
wave function is given by 



*(fc,P) = 



^ — P — 771 + ie 



(35) 



The valence wave function can be found by project- 
ing the Bethe-Salpeter wave function onto the light-cone 
x~^ = 0, see e.g. ]20| . Using the normalization convention 
of ElJ, we have 



7p{x, k,^i; A, A ) = — — / ^ 7+ 



2P 



2lT 



^ WA'((l-x)P+,P-L-k-' 



(36) 



where x is the fraction of the pion's plus momentum car- 
ried by the quark (x = k+ /P^), and the relative trans- 
verse momentum is k^^j = k^ — a;P^. The valence wave 



function is found from Eq. (M) to be 



V'(x,k^;A,A') ^ ■'^V^,^ Dw(x,k^ \Ml) 



x{\ — x) 



k-\ 5x\' - X777 (5a - A' 



(37) 



where we have employed the notation kx — + iXk^. 
As a result of the contour integration, we have a factor 
of 6[x{l — x)] implicitly in Eq. (p7|). Additionally the 
wave function is symmetric under interchange of x and 
1 — X. As is known, the simplistic form of this wave 
function leads to divergent quark distributions and form 
factors which will be handled below. Since this model is 
non-renormalizable, the choice of regularization scheme 
influences the dynamics. 



B. Form factor and generalized parton distribution 

The pion electromagnetic form factor for this model 
can be calculated from the Feynman triangle diagram. 
In order to extract the GPD, however, we need to choose 
the kinematics specified in Figure || with k as the mo- 
mentum of the struck quark. As it stands, using the 
wave function Eq. (|3^), the triangle diagram diverges. 
Following the approach of we covariantly smear the 
point-like photon vertex in Figure |l| in a way reminiscent 
of Pauli-Villars regularization 



^ A — 



A2 



A2 



fc2 - A2 + ie^' (fc + A)2 - A2 + ie' 



(38) 



Eq. (|3q ) is a simple way to model non-qq components 
of the wave function. Alternatively one could smear the 
qq — TT vertex in a covariant manner 1 23 , 2^ . This smear- 
ing should additionally respect the qq symmetry of the 
vertex. We do not pursue this option here since positiv- 



ity constraints (see section IV E) are likely to be violated. 
On the other hand, one could use Pauli-Villars subtrac- 
tions to regulate the theory, however, positivity would 
also be put into question. Because our concern is with 
model comparisons not phenomenology, we shall choose 
A = TO for simplicity. Although not obvious from in- 
spection, results for A 7^ to exhibit the same features 
investigated below. Most noteworthy, positivity remains 
satisfied when A 7^ to. 

Considering matrix elements of the current operator 
J'^ between pion states, the model ( p8[ ) conserves current. 
This can be demonstrated most easily by calculating A - J 
in the Breit frame. Additionally since the model is fully 
covariant, we can extract the electromagnetic form factor 
from any component of the current. In particular, poten- 
tial zero-mode singularities present in matrix elements of 
have been removed by the photon vertex smearing 
Eq.(^ 1^. Using the plus-component of the current, we 
have the expression 
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Tr 



ijfc + to) 75(|i- - p + m) + A + to) 7+ 



■[(fc + A)2 



ie] [(P-fc)2-TO2 + ie] 



(39) 



where the momentum transfer is t — A'^ and the skewness C is defined relative to the initial state: A+ = — C-P^ < 0. 
To calculate the GPD, we follow the procedure described in section |l^. The result can be written as 



(1 - C/2)H{X, C, t) = e{x - C)H,s{x, C, t) + 9{C - X) (ilinst(x, C, t) + ffnvai(a;, C, t)) • 



(40) 



where H^s is the piece determined by the effective two-body wave function, i?inst is the contribution from instantaneous 
propagation of the spectator quark, and the remaining contributions we term non- valence (although strictly speaking 
the instantaneous piece is also of the non- valence variety). It is a peculiarity of this model that explicit instantaneous 
terms are not present for x > 
The functional forms are 



7^-T3 X! ^osix', k'-^; A, A')V'off(a;, k^; A, A'), 
^ ""-^ A, A' 



4CDw'(x^k^|t) 
(1-x) x"(l-x") 

2(k^ • k'^ +to2) Dw(.T, k^ |M2) Dw'(x", k"^\t) 



x{l - x) x"(l - x")a;'(l - x')(l - C) 
X 2C Dw(x", k"^|t) + Dw(x, k^ |M2)] , 

where we have defined the effective wave function 



V'cff(2:,k-^; A, A') = 



g^/K/2Cm^ 



xHl 



A "A, A' 



Ato 5 



X- A' 



Dw'(x,k^|Af2) 



(41) 
(42) 

(43) 
(44) 



and made the abbreviation A = g'^\C\'^Ncrn^ /2{2tt)'^ . It is sensible to think of Eq. (^4|) as an effective wave function 
since x ^ 1 ~ x symmetry has been lost. Moreover the wave function vanishes at a; = 1 and is non-vanishing at 
a; = 0. This is a desirable addition to the dynamics stemming from the regularization. Notice the ladder kernel ( p^ is 
momentum independent and hence does not vanish at a; = 0, 1. This is the dynamical reason why the un- regularized 
wave function Eq. ( p7| ) does not vanish at the end points. Continuity of the GPD at a: = C follows directly from 
Eqs. dUil)- 

In Figure |lV^, we plot the GPD for the parameters: M = 0.14 MeV and to = 0.33 MeV. On the left, the graph 
shows the GPD for a few values of C as a function of a; for fixed while on the right we have fixed ^ and t varying. 



C. Sum rule and polynomiality 

Since not manifest, one should check the covariance of the model Eq. (^). With the covariant starting point 
Eq. (B9), we anticipate polynomiality will be satisfied which provides a useful check on our expressions Eqs. (ETl-Bs). 
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FIG. 4: Polynomiality conditions checked for the GPD 
Eq. @. The moments P„(C, t = 0) from Eq. @ are plotted 
as a function of for n = 0, 1, 2, 3. Additionally X's denote 
the simple (n + l)-point polynomial fit to the moment Pn- 
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FIG. 5: Comparison of covariant GPDs for the spinor triangle 
diagram. The GPDs Eq. @ (denoted LC) and Eq. (0) 
(DD-based) are plotted as a function of x for fixed = 0.9 
and t = — 4m^ for the mass = 0.15m. We also plot the 
difference between the two curves (5). The area under the 
curves is identically F{-4m'^) for LC and DD-based GPDs, 
and hence zero for their difference S. 



First we define the moments of the GPD with respect to 
asymmetric variables 



PniCt) - J x''H{x,C,t)dx. (45) 
Polynomiahty requires the moments Pn to be of the form 



(46) 



The zeroth moment is merely the sum rule for the form 
factor, hence aoo{t) — F{t). For simplicity, we check a 
few of the lowest moments for the polynomiahty condi- 
tion Eq. ( ^ ) at i = 0. In Figure [VB, we plot the mo- 
ments Pn{C, 0) for n = 0, 1, 2, 3 which appear as smooth 
functions. Additionally we plot simple (n-)-l)-point poly- 
nomial fits to the moments, which line up nicely with the 
integrals (p5|). 



D. Double distribution 

To construct the DD, we shall first proceed incorrectly 
by appealing to Lorentz invariance as in section l|. This 
will at least lead to one component of the DD, which can 
be compared to Eq. (pO|). 

Using Eq. (^) in the A+ = (Droll- Yan) frame, we 
can write F{t) in the form ([l^) with 

F{x, y; t) = {Zn? ~ MIx{1 - x) + y{\ - x - y)i) 
2nA 9{l-x~ y) y{l-x- y) 



(1-x) 



KQx(l -x)- y{l ~x- y)t 



(47) 



Aside from factors arising from spin, this DD is basically 
the same as that considered 



by utilizing = — M^/4 + m? . Not surprisingly, then, 
this DD satisfies Eqs. (pl]-p^. For reference we give the 
quark distribution function 



q{x) 



2^_4 (1 ~ xf 3m2 - Mlx{l - x) 



6 



- M^xil - x)Y 



(48) 



which one can realize 



which could be calculated directly from ijjcs in Eq. 

In Figure |, we plot the GPD Eq. (|^ as well as the 
GPD derived from DD Eq. ( [47| ) via Eq. (|l|). As in sec- 
tion U, the two are different despite the fact both models 
are covariant and posses the same form factor and quark 
distribution. Additionally we plot their difference as a 
function of x. 



E. Positivity Constraints 

Here we demonstrate another difference between the 
GPD in (^ ) and the one stemming from the the one- 
component DD Eq. (p7|). To do so, we look at the posi- 
tivity constraints. Originally these constraints appeared 
in jl^, ^ and were derived from the positivity of the den- 
sity matrix by restricting the final-state parton to have 
positive plus-momentum (and ignoring the contribution 
from E{x,(,t) for the spin-i case). Although the ma- 
trix elements involved for GPDs are off diagonal, they 
are still restricted by positivity and their diagonal ele- 
ments. Correcting the constraints for the presence of 
the i?-distribution was first done in [Q. By consider- 
ing the positivity of the norm on Hilbert space, stricter 
constraints for the spin-^ distributions H and E have re- 
cently appeared as well as constraints for the full set of 
twist- two GPDs Hi. 

For the scalar pion case there is of course no contri- 
bution from the non-existent i?-distribution and hence 
the original bounds are actually correct (modulo factors 
due to the difference of a spin-i proton versus a spin-0 
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FIG. 6: Comparison of GPDs: GPD calculated from the DD 
Eq. @ via Eq. ( p^ ) (DD-based) compared with the light- 
cone projection of the form factor Eq.(^) (Light-cone) for 
fixed C = 0.4 at t — 0. Here we plot the ratio R{x, (") appear- 
ing in Eq. (^^ as a function of x > (. Positivity constrains 
this ratio to be less than one. 



where q{x) is the model distribution function in Eq. (|4S^ 
Of course the above result holds for finite —t, however 
since the function T decreases with — t, Eq. ( pO| ) is the 
tightest constraint. Notice for ^ 0, the limit t = is 
in an unphysical region. If we treat this limit as formal, 
however, and analytically continue our expressions, we 
can use Eq. (|5C|). Such continuation is consistent with 
the light-cone Fock space representation of GPDs ^. 

Given the constraint Eq. (|50|), we can test whether 
GPDs calculated from the light-cone projection ( |4l| ) and 
DD ( ^ ) satisfy positivity. In Figure 6, we plot R{x,() 
for each GPD as a function of x for the fixed value oi ( = 
0.4. There is noticeably different behavior in the figure: 
positivity is violated by the DD-based model. As above 
(section [II ) , we must carefully derive contributions from 



the other component G{x,y;t). 



pion). Given the matrix element definition of the GPD 
consistent with equation (||), namely 



1 r dz- 
1-772 7 1^ 

x(7riP'mO)^+ijiz-)\niP)), (49) 



the spin-0 positivity constraint (for x > reads 

y'q{x}q{x') 



(50) 



Derivation of the correct DDs 



To derive both F and G DDs for the spin-i pion 
model, we must consider the action of the operator 

j\-^iD^'^ ■ ■ ■ iZ3^"] between non-diagonal pion states. In- 
serted into Eq. ( ^9[ ) which is now taken in the symmet- 
rical frame, this operator produces 



i-^ f 4 Tr{(^ + m)75(|^-;^-KA/2-l-m)75(|t + A + m)7^}(fc-f A/2)''i •••(fc-l- A/2)^"] 
^ J [fc2-m2-f ie]^[(fc-f A)2-m2 + ie]^[(fc-P-fA/2)2-m2+ie] 



(51) 



The presence of the trace 



Tr{. . .}^ = 4 Pf'im^ - P - k ■ A) + ^{m^ - + 2k ■ P) + k'^im^ - k^ ~ ^ + 2k ■ P ~k ■ A) 



(52) 



complicates evaluating Eq. (pl|) by requiring contributions from diagrams reduced by one propagator. Since we 
have smeared the photon via (pq) , the reduced diagrams are finite. Let us denote the propagators simply by 2t = 
(fc — P -|-A/2)2 — m'^ + ie, ^ — {k + A)2 - 777,2 -f ie ^^^^ <L — k"^ — m"^ + ie. To correctly evaluate Eq. (|5l|), we must 
write the trace as 



Tr{. . r = 4[f^(^ - i(<B + £)) + ^(Af2 - ^ - 21 + i(» - €)) + k^{M^ - ^ - ^l)^ 



(53) 



and evaluate each term separately canceling propagators in the denominator of Eq. (pl|). These integrals can easily 
be evaluated using Feynman parameters. For example, let us consider the non-reduced contribution from Eq. ([5^). 
The denominator of Eq. ^ appears as 21 S2 £2 g^j-^^j gQ introduce two Feynman parameters {x, y} to render the 
denominator specifically in the form [a;2l + j/*B + {1 — x — y)(t\~^ . One then translates k^^ to render the integral (hyper-) 
spherically symmetric via the definition fc^ = l'^ + /3P ~ {a+ 1) A'^/2. Here (3 = x and a = x + 2y—\. Using a Wick 
rotation to evaluate the resulting integral over /, we can cast the contribution to Eq. ( |5l| ) from non-reduced terms in 
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the form 



t:A I dp f '^da][{l-l3f-a^]D{f3,a;tf(2P{f3M^ + {l~p)t/2)-Aa{A4-t/2)^^'' 



-i+fi 



E 

fc=0 



7.1 



k\{n-k)\ 



(54) 



r 



with D{P,a;t) given by Eq. Given the form of 

Eq. (p4|), we can identify a} as DD variables and 
hence read off contributions to F and G DDs. 



SF{P,a;t) = t:A [(1 

X (/3M2 + (1 - f])t/2) 

5G{(3, a- t)^^Aa [(1 - [if - a"] D{p, a; 

ri . (56) 



Notice that these contributions respect the properties of 
DDs, namely bF is even in a while (5G is odd. This need 
not be the case, however, for each intermediate step of 
the calculation, e.g. contributions from Q3-reduced terms 
and ^-reduced terms are individually neither even nor 
odd in a while their sum is even and difference is odd. 

Ignoring for the moment contributions from 2l-reduced 
terms, we arrive at the DDs 



F{P,oi;t)^Ti AD{P,a;tf (1-/3)™' 
-/3a2Af2 + (1-/3) [(1-/3)2 



]i/4 



G(/3, oi;t)^~'K Aa £>(/?, a; tf 



(57) 



(58) 



The contribution from 2t-reduced terms has the form of a 
D-term. Using Feynman parameters for the denominator 
03^ £^ and suitable changes of variables, we arrive at the 
contribution to Eq. ( [sTI ) 



da- 



7r^a"+i(l-a2) 
«2 _ (1 _ a2)i/4]^ 



(-tr-(-f) 



A\M..] 



from which we can identify the _D-term 



(59) 



Switching now to asymmetric variables, we have 
F{x,y]t) ^ n A D{x,y;tf (1 - a;)m^ 

-x{x + 2y- If Ml + (1 - x)y{l - x - y)t 
G{x, y; t)^-TiA {x + 2y-l) D{x, y; tf 

+y{l -x-y)t- Ml{l ~x-{x + 2y- if) 
D{y; t) = nAy{l - y){2y - 1) [m^ - y{l - y)t] 



(61) 



(62) 



(63) 



where the function D{x,y;t) is given by Eq. (|29|). Ac- 
cordingly F{x, y; t) is Miinchen symmetric and G(a;, y; t) 
is antisymmetric, while D{y] t) is antisymmetric about 
y = 1/2. Notice F{x, y; t) in Eq. ^) is not that of 
Eq. (p^). Given the ambiguity inherent in defining F 
and G DDs (cf. Eq. ([2^)), there is no reason to believe 
the Fs would be the same. In principle, we could con- 
struct a gauge transformation [|5| to render the Fs the 
same. This would enable ident ification of the missing G 
function unique to section IV E . We shall not pursue this 
tangential point. ^ 

The function F(x, y; t) satisfies the reduction relations: 
it reduces to the quark distribution via Eq. (|l|) and 
integrates to the form factor (^3|) — the latter can only 
be checked numerically. Lastly then it remains to see 
whether the DD-b ased G PP lines up with true GPD cal- 
culated in sectio n |IVB . To construct the GPD we use 
the form of Eq. (|32|) modified to handle the Z?-term in 
Eq. (|63| ) separately 

H{x,C,t)^ dz dyS{x-z-Cy) 



[^"(^-^")li.(x";.) 



Fix,y;t) + ^^Gix,y;t) 



2(2-0 



(64) 



In Figure 0, we plot the GPD Eq. (|^ and the DD-based 



Although strictly speaking, the D-term is a contribution 
to the G-DD, we shall treat it separately for ease. 



^ Notice the contribution to the GPD from the D-term resembles 
that of i^inst in Eq. (kd) but is not identical. Both terms originate 
from a reduction of trie spectator's propagator. In the case of the 
D-term, the spectator's propagator is completely removed by the 
2l-rcduction. For i?inst, there is residual x-dependence stemming 
from the light-cone instantaneous propagator 7+ /2P+(1 — x). 
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FIG. 7: Contributions to the covariant GPD from DDs. The 
light-cone GPD Eq. and the DD-based Eq. @ are iden- 
tical, denoted (LC) and plotted as a function of x for fixed 
^ = 0.9 and t — —Am? for the mass M-^ = 0.15m. We also 
plot the individual contributions from F in Eq. (^ij) and from 
D in Eqs. ^ and denoted (G + D-term). 



2-C ^ 2(2-0 



Eq. (Q) . They are identical. We also plot the individual 
F and G + D contributions to the GPD. Even in the 
absence of the D-term, the contribution from G cannot 
be neglected in ascertaining the DD. 

The contribution to the GPD from F in the Figure 
is markedly different from Eq. ( ^o| ) and the analogous 
contribution from Eq. (^). These three GPDs, however, 
share the same form factor and quark distribution. Thus 
one can use DD-gauge freedom ||l5| to transform F into a 
new function and throw away contributions from G. The 
result is an infinite set of different GPDs with identical 
form factors and quark distributions. As pointed out in 



section IVE, one must be careful to maintain positivity 
although it is likely that there still is an infinite set of 
GPDs which would. 



V. SUMMARY 

Above we consider two covariant models for pions: one 
with scalar constituents and the other with spin-i. The 
spin-i model requires regularization and we choose the 
method in |^ in order to maintain positivity. For each 
case we derive the GPD from its matrix element defini- 
tion which forces us to consider the triangle diagram for 
the form factor with the plus momentum of the struck 
quark kept fixed in a general A''" 7^ frame. 

We also construct the DDs for each model. The ap- 
proach of ijll] leads only to one component of the DD (the 
"forward- visible" piece), thus resulting GPDs are incor- 
rect and need not satisfy positivity. This fact remains 
true even for C-odd distributions. The "gauge freedom" 
inherent in defining F vs. G DDs could be exploited, 
however, for phcnomcnological studies. To obtain the 
other component of the DD, we calculate the matrix ele- 
ments of twist-two operators. Resulting DD-based GPDs 
then agree with those calculated on the light cone. 
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